This is a preprint of a paper whose final and definite form will appear in Cent. Eur. J. Phys. 
Paper submitted 30-Jan-2013; revised 12-March-2013; accepted for publication 15-March-2013. 



This is a preprint to appear in Central European Journal of Physics. 



Noether's theorem for fractional variational problems 
of variable order 

Research Article 

Tatiana Odzijewicz^ *, Agnieszka B. Malinowska^ ' ^ ^ Delfim F. M. Torres-*^ -I- 

1 CIDMA — Center for Research and Development in Mathematics and Applications, 
Department of Mathematics, University of Aveiro, 3810-193 Aveiro, Portugal 

2 Faculty of Computer Science, Bialystok University of Technology, 15-351 Bialystok, Poland 



Abstract: We prove a necessary optimality condition of Euler-Lagrange type for fractional variational problems with 

derivatives of incommensurate variable order. This allows us to state a version of Noether's theorem 
without transformation of the independent (time) variable. Considered derivatives of variable order are 
defined in the sense of Caputo. 

PACS (2008): 02.30.Xx, 45.10.Db, 45.10.Hj 

Keywords: variable order fractional integrals • variable order fractional derivatives • fractional variational analysis 
• Euler— Lagrange equations • Noether's theorem 



1. Introduction 

Fractional calculus is the scientific discipline that deals with integrals and derivatives of arbitrary real (or complex) 
order. Since the XVIIth century, when fractional integration and differentiation was brought up for the first time, 
many well-known mathematicians contributed to the theory, among them Euler, Laplace, Fourier, Abel, Liouville 
and Riemann. For a comprehensive knowledge of fractional calculus we refer the reader to the books [20, 21, 30, 37]. 
The XXth century has shown that fractional order calculus is more adequate to describe real world problems than 
the integer/standard calculus [9, 18]. Therefore, not only mathematicians have currently a strong interest in the 
fractional calculus but also researchers in applied fields such as mechanics, physics, chemistry, biology, economics, 
control theory and signal processing [6, 29, 39]. 

A generalization of the fractional calculus was proposed in 1993 by Samko and Ross [38]. They considered integrals 
and derivatives of order a, where a is not a constant but a function. Afterwards, several works were dedicated to 
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these operators [35, 36]. Interesting applications of the proposed calculus are found in the theory of viscous fiows 
and meclianics [10, 12, 23, 31, 32]. 

In 1996, Riewe initiated the theory of the fractional calculus of variations by considering problems of mechanics 
with fractional order derivatives [33, 34]. Nowadays, the fractional calculus of variations is under current strong 
research (see, e.g., [1, 2, 5, 7, 8, 11, 15, 16, 25, 26] and references therein). For the state of the art on the subject 
we refer to the recent book [24]. Here we remark that results for problems of the calculus of variations with 
variable order fractional operators are scarce, reducing to those found in [4, 27, 28]. In particular, no Noether's 
symmetry theorem of variable order has been proved before. 

The main aim of the current work is to generalize the Noether theorem, originally proved by the German mathe- 
matician Emmy Noether in 1918, asserting that invariance properties of integral functionals lead to corresponding 
conservation laws [19, 22, 40]. In contrast with previous works [3, 13-15], where fractional versions of the Noether 
theorem for a constant non-integer order a are obtained, here we consider more general fractional variational 
problems with variable order derivatives. 

The paper is organized as follows. In Section 2, basic definitions and properties of variable order fractional 
operators are given. Sections 3, 4 and 5 contain our main results: we prove integration by parts formulas for 
variable order fractional integrals (Theorem 3.1) and derivatives (Theorem 3.2), a necessary optimality condition 
of Euler-Lagrange type for fractional variational problems with derivatives of incommensurate variable order 
(Theorem 4.1), and a version of Noether's theorem for such fractional variable order problems (Theorem 5.1). 
Then, in Section 6, we illustrate our results through an example. We finish with Section 7 of conclusion. 

2. Variable order fractional operators 

In this section we recall the basic definitions necessary in the sequel, and derive a new interesting mapping property 
for the fractional integral operators of variable order (Theorem 2.1). We set A ~ {{t, r) G : a < r < t < 6}. 

Definition 2.1 (Left and right Riemann— Liouville integrals of variable order). 

Let a : A ^ (0, 1) and / G Li [a, h\. Then, 



t 




t > a, 



a 



is called the left Riemann-Liouville integral of variable fractional order a{-,-), while 



b 




t < b, 



denotes the right Riemann-Liouville integral of variable fractional order a(-, •). 
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Theorem 2.1. 

Let a{t, t) = a{t — r) with < a{t — r) < 1 — — for all < t — r < b and a certain n £ N greater or equal than 
two. If f £ AC[Q,h], then 



ol]~-^-'-^ fit) 



r(l-Q(t-r)) 



(t-r)-"('-"'/(^)rf^ e AC[Q,b]. 



If f £ AC[-b,Q], then 



U 

t 

Proof. We give the proof for the left integral; the other case being proved similarly. Let k{s) 
Since < a(s) < 1 - 

1. for s < 1 we have Ins > and s""'"' < 1, 

2. for s < 1 we have Ins < and s""'"' < S"~^. 
Therefore, 



r(l-a(s))' 



-a(s) 



/ lfc(s)|dr: 
JO 



-ci(s) 



r(i-Q(s)) 



dr < 



r(l-a(s)) 



ds 



Using the inequality 



x"^ + 1 



valid for all x £ [0, 1] (see [17]), we obtain that 



(1) 



-ds + 



r(l-a(s)) A r(l -a(s)) 



— s " ds < I ds + s " ds = 1 + " — n < oo. 



J 1 



It means that k{s) G ii[0, 6]. Moreover, the condition / G AC ([0, fe]) implies that 



fit) = / + /(O), where g e Li ([0, 6]) . 

Jo 



Define 



and integrate: 



hit) := 



r(l-a(s)) 



s "^"^ git — s)ds 



fhie)de= f de f ^ .-'^(^^gie - s)ds. 

Jo Jo Jo r(l"a(s)) 



(2) 



3 
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Observe that 



r(i-Q(s)) 



d9 ds = 



-a{a) 



<\\gh 



< oo. 



r(l-a(s)) 

b 1 



1^ \g{9-s)\de^^ds 



r(l-a(s)) 



6-"("' ds 



It means that assumptions of Fubini's theorem are satisfied, so we can change an order of integration in (2). 
Hence, 

I' Hoyo ^ I' ds [ ^(T^-- - s)de ^ I - s)de. 

Putting ^ = 6 — s, one has = d9 and we get 



Jo 



h{e)de = 



/o r(i-Q(s)) 

But " giO'^^ ~ /(^ ^ ^ /(O) and, therefore, the following equality holds: 



h(e)de^ /*— ^——s-"^"^ fit - s)ds ~ f(0) r — L__s— ("jrfs. 



q(s)) 



Putting T = t — s, dr = —ds and 

io io r(l-«(t-r)) Jo r(l -a(t-r)) 

Both functions on the right-hand side of the equality 

Trn — 77 — ^(t - r)-"('-^'/(r)dr = fh{e)de + f{Q) f ^ -Jt-r) 

Jo r(l -a(t-r)) Jo Jo r(l -a(f-r)) 

belong to AC{[0,b]), hence oi^''^ '/ £ ^C[0,&]. 



dr. 



a(t-T) 



dr 



□ 



Definition 2.2 (Left and right Riemann Liouville derivatives of variable order). 

Let a : A — ^ (0, 1). If alt~'^^ ' f £ AC[a, b], then the left Riemann-Liouville derivative of variable fractional order 
a(-, •) is defined by 

a 

while the right Riemann-Liouville derivative of variable order a{-, •) is defined for functions / such that til g 
AC[a,b] by 



< b. 
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Definition 2.3 (Left and right Caputo derivatives of variable fractional order). 

Let < a{t, r) < 1 for all t,T £ [a, b]. If / € AC[a, b], then the left Caputo derivative of variable fractional order 
q(-, •) is defined by 

t 

a 

while the right Caputo derivative of variable fractional order a(-, •) is given by 



3. Variable order fractional integration by parts 

In this section we derive integration by parts formulas, which are essential for proving Euler-Lagrange equations. 

Theorem 3.1 (Integration by parts for variable order fractional integrals). 

// i < a(t, r) < 1 for all (i,T) G A and a certain n G N greater or equal than two, and f,g £ C ([a,fe];R), then 

f g{t).it^ '-^f{t)dt = f f{t),t ''^g{t)dt. 

J a J a 

Proof. Define 

F{r,t):-^ 



]p(^(t~rr(*--)-iff(t)/(r)] ifr<f, 
if r > 



for all (r, f) G [a,b] x [a,b]. Since / and g are continuous functions on [a,b], they are bounded on [a,b], i.e., there 
exist Ci,C2 > such that \g{t)\ < Ci and \ f{t)\ < C2, t e [a,b]. Therefore, 



dr dt 



Because ^ < a{t,T) < 1, 

1. for 1 < t - r we have \n{t - r) > and {t ~ r)°'*'^'"^ < 1; 

2. for 1 > t - r we have ln(t - r) < and {t - r)"'*'^)-^ < {t - r)""\ 
Therefore, 



C1C2 



< C1C2 



Nocthcr's theorem for fractional variational problems of variable order 



Moreover, by inequality (1), valid for x £ [0, 1], one has 



C1C2 



< CiC: 



CiC2{b-a) ( ^-l + n-a 



< 00. 



Hence, one can use the Fubini theorem to change the order of integration: 



J a 



□ 



Theorem 3.2 (Integration by parts for variable order fractional derivatives). 

Let < a{t,T) < 1 — ^ for all {t,T) G A and a certain n £ N greater or equal than two. If f £ ([a,6];R), 
g £ C {[a,h];W), and til "'''5, then 

r g{t)^Df-y{t)dt^ f{t)tll-''^-^git)\' + [' f{t)tDf-^g{t)dt. 

If f {[a,b];'R), g e C {[a,b]-R), and all'"^'''^ g e AC[a,b], then 

r g{t)?D^^-yit)dt^- f{tUl-"^---^git)\' + f f{t)aDf - ~>g{t)dt. 

Proof. By Definition 2.3, it follows that ^D"^'''^f{t) = a4 ' ' ^/(t). Applying Theorem 3.f and integration 
by parts for classical (integer order) derivatives, we obtain 

1^ g{tf^Df---\f{t)dt = I' g{t)J--^-'-^^J{t)dt = ±f^t)tfr^-'-^g{t)dt 

= mtfr^-'-'git)\^ - [ f{t)f^tfr^-'-^g{t)dt 

= f{t)ttr^-'-^g{t)\\ f f{t)tDf--^g{t)dt. 



The second formula is proved in a similar way. □ 
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4. Euler— Lagrange equations for incommensurate fractional 
variational problems of variable order 

Now we prove necessary optimality conditions of Euler-Lagrange type for incommensurate order fractional vari- 
ational problems. Let ai{t, t), Pi{t, t), i = 1, . . . , n, n G N, satisfy the assumptions of Theorem 3.2. Consider the 
following problem. 

Problem 4.1. 

Find a function q — q{t) for which the functional 

b 

a 

attains an extremum on the set 

V=[qe C\[a,b];R) : q{a) = qa,q{h) = q, and :r 0^'*'"'?, f G C([a, &]; R), z = 1,.. .,n} . 

We assume that L £ ([a, b] x R^""*"^; R); f m> di+2L is continuous, has absolutely continuous integral tl^~"''''^' 
and continuous derivative tD|^''''^' for each i = 1, . . . , n; f M- d„+i+2L is continuous, has absolutely continuous 
integral ^Jf and continuous derivative aOf''*'^' for each i — l,...,n. For simplicity of notation, we 
introduce the following notation: 

{q,a,l3}{t) (t,gW,fDr'*'^'g(t),fDf<*'^'g(t)) , 

where 

Crya(t,T) _ (Cjjai(t,T) Cjja„(t,T)\ C p,P(t,T) _ (Cp,l3l(t,T) C t-,/3„ (t .r) \ 

at ' [at ^ ' ' ' ^ a t /' t Jj • It fj j ' ' ' j t fj j ' 

Definition 4.1. 

Let / G AC{[a,b]). For til '''"'^'ff £ ^Cilo-jb]) we define the following operator: 

DI*''^' [/, 9] -m^'^^'g + 5a D?''-^' [/], (4) 
while for aft^''^''^^ g G AC{[a,b]) we define 

Df-^\f,g] := -/.-DZ'^^'M +3fD^'*-'[/]. (5) 

Remark 4.1. 

If 7(t, t) = 7 is a constant function, then 

hm Dl [/, g] = fg' + gf = ^(fg) = lim Dl [g, f] 

and 

lim Dl[f,g]^-fg' -gf ^-^ifg)= lim Dj^,/]. 
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Theorem 4.1. 

Let function q be a solution to Problem 4^.1. Then, 

n n 

&L {q, a, /?} (t) - ^ [1, d.+2L {g, a, /?} (t)] - ^ D'^''*'^' [1, a„+2+.L {g, a, /?} (t)] = 0. (6) 

i = l 1=1 

Proof. Assume that q is an extremizer of J . Consider the value of J' at nearby function q{t) — q(t) + erjit), 
where e G R is a smaU parameter and rj G C^([a,6];R) is an arbitrary function satisfying 77(a) = 77(6) = and 
such that ^D°''*'^^77 and pD^*'''''^'?) are continuous. Let 

J{e) = J[^]= f L{q,a,p}{t)dt. 
A necessary condition for q to be an extremizer is given by 

^ = ^ /' (d2L {q, a, /?} {t)ri{t) + ft+aL {q, a, /?} (t)^ D^^^''^^r,{t) 

+ a„+2+,L {g, a, /?} D^(*'^),7(t) j dt = 0. (7) 
Using the variable order fractional integration by parts formulas (Theorem 3.2), we obtain that 



b 



d.+2L{q,a,P}{t)^D:^^''^^r^{t)dt 



b 



v{t)tll-'''^'-'^^d,+2L{q,a,l3}{t)\ + / 77(t)tD°'(*'^'a+2L{g,a,/3}(t)df, t^l,...,n, 



and 



J a 



+i+2L{q,a,l3}(t)?D^,^^''^^rj{t)dt 



b 



= -r,(t),4-^'^''"'9„+,+2L{g,a,/?}(t)f + / v{t)aD^'^''^^ ^r.+^+2L {q, a, p} {t)dt, i = 1, 



, n. 



Because Tj{a) — Tj{b) — 0, (7) simplifies to 

/[, / n n \ 

V{t) d2L {q, a, /?} (t) + J2 tD"''''^'ft+2L {g, a, /?} (t) + ^ .Df''*'"^a„+.+2L {g, a, /?} (t) = 0. 
V i=i i=i / 

Applying the fundamental lemma of the calculus of variations (see, e.g., [40]), we obtain that 

n n 

d2L {g, a, /?} (t) + tD''/'-^'>d,+2L {g, a, /?} (t) + „Df' (*'^)a„+,+2L {g, a, /?} (t) = 0. 
Finally, by Definition 4.1, we arrive to (6). □ 
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5. Noether's theorem of fractional variable order 



Now we show, employing operators (4) and (5), that invariance of (3) conducts to a variable order version of a 
fractional conservation law. To prove this, we borrow the method from [13, 15, 19], where Noether's theorem is 
stated for fractional derivatives of constant order. 

Definition 5.1. 

A function q that is a solution to (6) is said to be a variable order fractional extremal for functional . 
Definition 5.2. 

We say that functional (3) is invariant under an e-parameter group of infinitesimal transformations 



q{t) = q{t)+ei{t,q{t)) + o{e) 



(8) 



L {t, qit),^Df'-'-^qit), f dt = j\{t, q{t),aDf'--^m, f "^'^(t)) dt (9) 

for any subinterval [ta,tb] CI [a,b]. 

Lemma 5.1 (Necessary condition of invariance). 

If functional (3) is invariant under an e-parameter group of infinitesimal transformations (8), then 

n 

d^L {q, a, P} {t)at, q{t)) + J2 d.+2L {q, a, p} (f)^Dr q{t)) 

n 

+ J2 ^"+'+2^ il' /^} q{t)) = 0. (10) 

1=1 

Proof. Since, by hypothesis, condition (9) is satisfied for any subinterval [fa, it] Q [o-,b], we have 



(t,q{t),^Df''-^q{t),?lf,^''-k{t)) = L [t,qit),aDf''-k{t),?lf,^''-^qit) 



(11) 



Now, differentiating (11) with respect to e, then putting e = 0, and applying definitions and properties of variable 
order Caputo fractional derivatives, we obtain that 



= d2L{q,a,P}{t)at,q{t)) 

n 

+ J29^+2L{q,a,p}{t) 



de 



+ ^a„+,+2L{g,a,/3} {t) 



„ r(l- a.(t,r)) 

b 1 



{t-ry 



dr 



q{T)dT 



de 



r(l-ft(t,r)) 



dr 



q{T)dT 



d2L {q, Q, /?} mit, q{t)) + J2 9^+2L {q, a, /?} (t)^Dr<''"'C(t, <?(*)) 

n 

+ dn+^+2L {q, a, p} (t)f ('■^'C(t, q{t)). 



□ 
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Theorem 5.1 (Noether's theorem for variable order fractional variational problems). 

If functional (3) is invariant in the sense of Definition 5.2, then 

n n 

^ D^''*'^' m q{t)),d.+2L {q, a, /3} (t)] + ^ D^'<''"'[e(t, q{t)),d„+.+2L {q, a, /?} (i)] = 0, t G [a, b], 

i=l i = l 

along all variable order fractional extremals q(-). 
Proof. By Theorem 4.1 we have 

n n 

d2L {q, a, P} {t) = -J2 tD"'^''^'9.+2L {q, a, /?} (t) - ^ „Df-'*-^'a„+.+2i {q, a, /?} (t). (12) 
Substituting (12) into (10), we obtain 

n n 

- ^(*' g(t))t-D^'^''"'a,+2L {g, a, /?} (t) - ^ g(t))aDf'(*'^'a„+,+2L {g, a, /?} (t) 

i=l i^l 

n n 

+ J2 d^+^L {q, a, P} {tf^D'i^ ('-^'^(t, q{t)) + ^ d^+^+^L {q, a, /?} (t)f ■ ?(*)) = 0. 

i=l i = l 

Finally, by Definition 4.1, one has 

^ D^'<*'^' K(t, g(t)), a+2L {g, «, /?} (*)] + ^ D^"'*'"' K(t, g(t)), a„+.+2i {g, a, /?} (t)] = 0. 

1=1 i = l 

□ 

6. An illustrative example 

Let a and /3 be two functions such that < a{t, t), I3{t, r) < 1 — j for a certain natural number I greater or equal 
than two. Consider the following problem: 

J[q] = f L (t, ^ f '*'"'q(t)) dt extremize 

J- ^ ' (13) 

g(a) = qa , q(b) = q^. 

For the transformation 

g(t) = q(t)+£c + o(e), (14) 

where c is a constant, we have 

L (t, f Dr'*-'g(t), f D«'-'g(t)) dt = L (t, f Dr'*-',-(t), f (*-',-(t)) 
for any [ta,tb] CI [a,b]. Therefore, J[q\ is invariant under (14) and Noether's theorem (Theorem 5.1) asserts that 

D^!(*-> [c,9.L(t,fDr(*-'g(t),fDf(*-'g(t))] + of [c,d,L {t,'=^Df'-~> q{t),<^ li^'-^~> q{t))\ = Q 
along any extremal q{t) of (13). 
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7. Conclusion 

Noether's symmetry theorem, establishing that variational invariance imply conservation laws, is one of the most 
important results in physics and the calculus of variations. Such conservation laws are, however, only valid for 
conservative systems. Nonconservative forces, like friction, remove energy from the systems and, as a consequence, 
Noether's conservation laws cease to be valid. In order to cope with dissipative forces that do not store energy, one 
possibility is to use fractional calculus [33, 34]. The study of problems of the calculus of variations with fractional 
derivatives of variable order is, however, a rather recent subject, and available results reduce to those available 
at [4, 27, 28]. Here we obtain Euler-Lagrange optimality conditions for variational problems with fractional 
derivatives of incommensurate variable order and, for such fractional Euler-Lagrange extremals, we provide a 
notion of conservation law and prove a version of Noether's theorem. 

The variable order calculus of variations is underdeveloped and much remains to be done. We mention here one 
important question that remains open. The Noether type theorem we obtain here only includes the terms related 
to "conservation of momentum". It would be important to obtain a more general form of the Noether theorem 
with "energy" terms. 
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